Color Dynamics In Phase Space: The Balescu-Lenard-Vlasov Approach 

O' 

O ■ A. Bonasera * "" 

O 

(N 



Oh: 

in 



o 



p. 



^ Laboratorio Nazionale del Sud, Istituto Nazionale di Fisica Nucleare, Via S. Sofia 44, 
95123 Catania, Italy . 



We propose to use the Balescu-Lenard-Vlasov (BLV) equation to describe relativistic 
heavy ion collisions. We use an inter-quark Richardson's potential consistent with the 
, ^ ■ indications of Lattice QCD calculations. The color degrees of freedom are explicitly taken 
i/-^ ! into account. We explicitly demonstrate that the Vlasov approach alone is insufficient 
O I in the hadronization region. In order to overcome this problem we prepare the initial 
Q I condition for many events using molecular dynamics with frictional cooling and a Thomas- 
O ■ Fermi approximation to the Fermi motion. These events are averaged and propagated in 
time using the Vlasov approach. The fluctuations are used to evaluate the collision term 
and in turn the number of q-q created. 



Oh! 

O , 1. INTRODUCTION 



One of the open problems in theoretical nuclear and particle physics is how to obtain 
k> ! the well known nuclear properties starting from the quark degrees of freedom [|l|. This 
H I also includes the possibility of understanding the basic free nucleon-nucleon interaction 
from quark and gluon dynamics. Some kind of solution to this problem is becoming more 
and more needed with the new experiments done or planned using ultra-relativistic heavy 
ions (rhic) at CERN and at RHIC The search for a quark-gluon plasma (QGP) in such 
collisions is in fact one of the new and most exciting directions in physics at the border 
between nuclear and particle physics 0. Quantum ChromoDynamics (QCD) because 
of its difficulties (numerical and conceptual), has been applied so far to some limited 
cases such as quark matter at zero baryon (p^) density and high temperature (T) ^^■ 
Furthermore in rhic dynamics plays surely an important role and accordingly the theory 
should be dynamical. 

Recently |@ , we have proposed a dynamical approach based on the Vlasov equation m| 
to reproduce hadron masses and the properties of nuclear matter at finite pB- Some works 
in the same spirit are discussed in 0,0. Our approach needs as inputs the interaction 
potential among quarks, which was borrowed from phenomenology i.e. the Richardson's 
potential]^, and the quark masses which were fitted to reproduce known meson masses 
such as the vr, the 0, the rjc etc. When the particles are embedded in a dense medium 
such as in nuclear matter (NM) the potential becomes screened in a similar fashion as 
ions and electrons in condensed matter do, i.e. Debye screening (DS)P,H. 
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Here we refine that approach in two important aspects which are the treatment of the 
color degrees of freedom and the inclusion of the collision term. In the previous works [Q 
color degrees of freedom were implicitly taken into account through the use of a Debye 
radius that effectively screens the qq interaction potential. In the present paper we give 
to the quarks explicitly a color (using the Gell-Mann matrices) and follow their dynamics 
in phase space solving the Vlasov equation (VE). Thus screening will be dynamically 
obtained. In general the self screening obtained in Vlasov dynamics is inadequate, which 
is the reason why it was not adopted in the earlier attempts 0. In fact we will show 
explicitly that the Vlasov approach alone gives a good description of the system at large 
densities only,i.e. in the QGP region. In order to overcome this problem we adopted 
the following strategy. We first prepare the initial conditions using molecular dynamics 
(MD) with frictional cooling for many events. The events are averaged and care is taken 
of antisymmetrization. These are the initial conditions for the Vlasov evolution. Since 
the VE fulfills the Liouville theorem, the initial phase space density remains constant in 
time. Thus the initial antisymmetrization and eventual clustering obtained in the cooling 
process are maintained. When simulating a nucleus-nucleus collision, inter quark corre- 
lations and particle production become very important. We have included this feature 
using the most general form of the collision term which is the Balescu-Lenard type. This 
is very crucial since we are dealing with particles which are interacting through long range 
forces. 

2. The Balescu Lenard Vlasov equation 

We outline our approach on purely classical grounds, however the same results can be 
obtained within the Wigner transform formalism [^,H of the quantum BBGKY-hierarchy 
in the limit h ^ 0. 

The exact (classical) one-body distribution function /i (r, p, t) satisfies the equation 
(BBGKY hierarchy) §: 

dth + ^ • VJi = j d{2)VrV{v, r^) Vp/2(r, v^, p, p2, t) (1) 



E = Jp^ + ml is the energy and rrii = lOMeV is the (u,d) quark mass. Here we assume 
the potential to be dependent on the relative coordinates only. A generalization to include 
a momentum dependent part is straightforward. /2 is the two-body distribution function, 
which in the classical limit reads: 

/2(r, r2, p, P2, t) = S^_^^5(r - r„)(5(p - p^,) x 5{r2 - r^)5(p2 - P/?) (2) 

where Q = q + q is the total number of quarks and anti quarks . Inserting this equation 
into Eq.(l) gives: 

dji + ^ • V./i - VrU ■ Vp/i = (3) 

Where U = T,jV{r,rj) is the exact potential. Let us now define /i and U as sums of 
an ensemble averaged quantity plus the deviation from this average: 

fi = fi + 6h;U = U + 6U (4) 



Substituting these equations in Eq.(3) and ensemble averaging gives: 

dth + ^ ■ Vrh - V,f/ ■ Vp/i =< V,<5f/Vp5/i > (5) 

where one recognizes in the Ihs the Vlasov term and in the rhs the Balescu-Lennard 
colhsion term P,p,pi]|. The mean-field is given by: 
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f/(r) = ^SevSjV(r,rj) (6) 

In agreement to LQCD calculations [|l],|I| the interacting potential V{r) for quarks is {Ti = 
1): 



and ID 

f{t) = l-A I ^—-——— (8) 

^ ^ J q [ln{q^ - 1)]^ + tt^ ^ ^ 

We fix the number of fiavors nj = 2 and the parameter A = 0.25GeV . In Eq.(7) we 
have added to the Richardson's potential the chromomagnetic term (ct), very important 
to reproduce the masses of the hadrons in vacuum. Since in this work we will be dealing 
with finite nuclei, the ct can be neglected, we only notice that with the parameters choice 
discussed here, the hadron masses can be reproduced by suitably tuning the ct term 0. 

The A" are the Gell-Mann matrices. From lattice calculations we expect that there is 
no color transport for distances of the order of 0.2 — 0.3/m, which are distances much 
shorter than the ones we will be dealing with in this paper. Thus we will use the Aa^g only 
commuting diagonal Gell-Mann matrices (Abelian approximation)]^. 

Numerically the BLV equation(5) is solved by writing the one body distribution function 
as: 

1 ^ 
h{r,p,t) = —Y^a,{t)5{r-n{tmp~pm) (9) 

ntp j 

N = Qritp is the number of such terms. Actually, N is much larger than the total quark 
number Q, so that we can say that each quark is represented by rifp = N^v terms called 
test particles(tp). Inserting Eq.(9) in the BLV equation gives the Hamilton equations of 
motion (eom) for the tp [|]. The weight factor fulfills the equation dtaiit) = V n^UV ki^ fi 



T\\ . Initially the weight factor is equal to 1 for the tp. When for a particular tp the 
weight factor becomes larger than 2 a new tp (quark) is created in the same phase space 
region. In the same way if the weight becomes less than zero for a given tp another tp with 
negative sign ( an anti-quark) is created. Notice that on average over time and ensembles 
the number of newly created quarks is equal to the number of anti quarks. Particular care 
is taken in order to have specific color conservation as well. We would like to stress now 
that in our picture the mean field and the collision term are calculated in a consistent 



way starting from a two body potential. This is indeed a quite general method which can 
be applied to other physical systems. 

Initially, we distribute randomly the tp in a sphere of radius R = tqbA^^^ (the radius of 
the nucleus) in coordinate space and pj in momentum space, tqb = (4:7—)^''^, A = Q/3 



and pb is the baryon density, pf is the Fermi momentum estimated in a simple Fermi gas 
model by imposing that a cell in phase space of size h = 2n can accommodate at most 
two identical quarks of different spins, flavors and colors. A simple estimate gives the 
following relation between the quark density 



Uq, and the Fermi momentum: 



n„ 



9q 3 

g^2 / 



(10) 



The degeneracy number Qq = nj x UcXUs, where Uc is the number of colors and n^ is the 
number of spins [§. For quarks and anti quarks 3 different colors are used red, green and 
blue (r,g,b) i. 

3. Results 
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Figure 1. Energy per nucleon (top) and energy density of the quarks (bottom) vs. 
baryon density. 



In figure(l), we plot the total energy per nucleon (top) and energy density (in units of 
the Fermi gas energy density]^) vs. baryon density. The full triangles give the results 
obtained by randomly distributing the tp as described above. We notice that a minimum 
at about pc = 2.08fm~^ is found with E^/A = 2GeV/A. Such a minimum is at much 



higher density and energy than expected for the ground state (gs) of a nucleus (po = 
0.16/m~^ and Et/A = 0.938 — 0.016GeV/A). An important property of the system that 
we have described above is the following. If we rotate the quarks in color space, regardless 
of their position in r-space, the total energy will remain the same. This is indeed a "pure" 
Vlasov solution and demonstrates explicitly the in-capability of the Vlasov approach to 
give clustering of quarks into nucleons. However this result is already instructive since 
it gives us an hint on where the quark and nuclear matter are located, i.e. above and 
below pc respectively. This result is qualitatively in agreement with the Hartree-Fock(HF) 
calculations of refs. |[T^ , P^ (compare to fig.l in |jl2[)- For a discussion on why the Et/A 



increases at low densities in the HF/VE approaches we also refer to [|T2 . 

Of course distributing randomly the quarks in a sphere in r and p-space is not the most 
economical way to find the real gs of the system. In MD one searches for a minimum energy 
by introducing a friction term. The friction acts in such a way to lead the particles to a 
configuration for which the potential energy is a minimum. We cannot use this technique 
for our system since we are dealing with fermions and the friction term will destroy the 
initial antisymmetrization. In order to overcome this difficulty we adopt the following 
strategy. First we prepare iVe^ events as above, and we evolve them numerically solving 
the eom but with friction included. Because of the large number of particles interacting 
with attractive and repulsive forces, the quarks will slowly evolve to new positions where 
the potential is lower while keeping the initial root mean square radius approximately 
constant. When the averaged potential (over events) reaches a given value Vmin, we look 
for the two closest particles (j, k) to a quark i in the same event. For these three quarks 
we know what the local density and the number of colors are. For instance if in a certain 
region we find two red and one blue quark, we use n^ = 2 in Eq.(lO) and calculate the 
local density from the knowledge of the distances of the 3 quarks. In this way the Fermi 
momentum is defined locally similar to the procedure used in nuclear or atomic physics 
(Thomas- Fermi approximation) [|14|. We repeat this for all quarks i in all events and 
calculate the total energy for this state . We let the system evolve again with friction 
included to a lower potential energy Vmin = Vmin — ^V, where 6V is a constant. We 
calculate the local density and local color numbers again and apply the Thomas-Fermi 
approximation to obtain the new total energy. We stop the procedure when the total 
energy is a minimum. The initial conditions so obtained are then propagated in time 
using the VE in order to maintain the initial antisymmetrization. 

The open triangles in fig.(l) are the result of the minimization. Notice especially at 
low densities the large decrease of the total energy of the system as compared to the 
Vlasov result. Now the calculated total energy at the nuclear gs is very close to the 
experimental value indicated by the full circle. However, we find slightly lower energies 
at lower densities, i.e. the gs of the nucleus is shifted in our calculations at about 1/10 
of the experimental one. This should not be surprising in view of the simple potential 
that we have used. Also we have not tried to look for a best set of fitting parameters in 
these exploratory studies. At low pB, the global invariance for rotations in color space is 
lost, i.e. if we exchange the colors of two distant quarks, the total energy of the system 
will change. At high densities (larger than 2fm~^) the Vlasov and MD solutions are 
the same. This can be also seen in the bottom part of the figure where energy densities 



are given. We would like to stress again the qualitative agreement to ref.|T^ where a 



stochastic method had been used to calculate the g.s. energy of the system. This is quite 
evident if one compares our fig.l to fig.l of ||12|| . 

The energy density displayed in fig.l (bottom) is a smooth function apart some fluctua- 
tions around 2fm~^ density. From this result we can exclude a first order phase transition 
but a second order phase transition might be possible, see the discussion in |l^ . 
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Figure 2. Ne-I-Ne collision at 100 GeV/A and b=Ofm. The quarks are given by the 
(r,g,b) while the created anti quarks by the yellow, magenta and cyan colors. 



Once we have devised a method to find a (approximate) nucleus ground state starting 
from quark degrees of freedom we can perform nucleus-nucleus collision simulations. In 
figure (2) we plot snapshots of a simulation for Ne-I-Ne collision at 100 GeV/A. Initially the 



quarks are distributed randomly in a sphere of radius R and the minimization procedure 
is apphed. Such an initial condition is Lorentz boosted and the dynamics is followed in 
time. In this simulation N^.^ = 9 for numerical reasons and one event only is displayed 
in the figure (2). Because of the large fluctuations from event to event, the collision term 
is very important and new quarks and anti quarks are created. These are displayed in 
the figure with different colors. The newly created quarks are strongly sensitive to the 
parameters entering the two body force and especially to the linear term. This might have 
some connection with the Schwinger particle production mechanism |^ but more studies 
are needed to find a precise link (if any) . 

4. Summary 

In conclusion in this work we have discussed microscopic Vlasov/MD approaches to 
finite nuclei starting from quark degrees of freedom with colors. In order to obtain the 
correct initial conditions we have introduced a method based on MD with frictional cooling 
plus a Thomas-Fermi approximation for the Fermi motion. We have shown that the 
method is able to describe at least qualitatively the well known features of nuclei near 
the ground state. At high densities a second order phase transition from nuclear to quark 
matter is predicted. Such a transition is due to the restoration of global color invariance 
at high densities and we have defined an order parameter accordingly|T^. We have used 
the initial conditions for the nucleus ground state to simulate heavy ion collisions at 
ultra-relativistic energies with the Balescu Lenard collision term. The mean field and the 
collision term are calculated starting from the same elementary interactions. Collisions 
are connected to deviations from the average mean field. This gives also large fluctuations 
in the distribution function which results in new particles creation. Our approach can be 
very useful for the understanding of the quark gluon plasma formation and its signatures. 

I thank Toshiki Maruyama for useful discussions, suggestions and the use of drawing 
subroutines. 
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